Finite volume effects and dynamical chiral symmetry breaking in QED3 
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We investigate the impact of finite volume effects on the critical number of flavours, Nf, for chiral 
symmetry restoration in QEDg. To this end we solve a set of coupled Dyson-Schwinger equations 
on a torus. For order parameters such as the anomalous dimension of the fermion wave function 
or the chiral condensate we find substantial evidence for a large dependence on the volume. We 
observe a shift in Nf from values in the range of 3.61 < A'^^ < 3.84 in the infinite volume/continuum 
limit down to values below Nf < 1.5 at finite volumes in agreement with earlier results of Gusynin 
and Reenders in a simpler truncation scheme. These findings explain discrepancies in Nf between 
continuum and lattice studies. 
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I. INTRODUCTION 

Dynamical chiral symmetry breaking in the theory of 

Saantum Electrodynamics in (2+1) dimensions, QED3, 
has been studied quite intensively over the years. The 
problem is of considerable interest for two reasons. On 
the one hand, QED3 has enough similarities to Quantum 
Chromodynamics (QCD) such that analogies to the more 
complicated non-Abelian case may be drawn. On the 
other hand, QED3 itself is of considerable interest due to 
possible applications in condensed matter systems. 

In particular, QED3 with Nf — 2 flavours of mass- 
less fermions has been suggested as an effective low en- 
ergy theory of high- Tc cuprate superconductors [3, 0] ■ 
These possess an unconventional d-wave symmetry of the 
pairing condensate, with nodes at the electronic Fermi 
surface allowing for a description in terms of massless 
nodal quasi-particles. The quasi-particles do not couple 
to an external electromagnetic field and hence represent 
pure spin degrees of freedom; a spin-charge separation 
has taken place. The interaction of these spinous with 
the collective, topological excitations of the gap can be 
described by a U{1) gauge theory. Since furthermore the 
motion of the quasi-particles is mainly confined to the 
two-dimensional copper-oxygen planes in these systems 
one ends up with quantum electrodynamics in (2+1) di- 
mensions. Depending on whether the system is ordered 
or disordered it is either in an insulating quantum antifer- 
romagnetic (AF) phase or a pseudogap (PG) phase with 
remnant properties of the underlying superconductor. In 
the gauge theory the AF-phase corresponds to a phase 
with broken chiral symmetry and long-range correlations 
due to massless photons. In the PG-phase the fermions 
are massless and the fermion wave function as well as the 
photon propagator develop power laws at small momenta 
with a fractional anomalous dimension ^, ^] . 

The above considerations explain the interest in deter- 
mining Nj, the critical number of fermion flavours for 
the chiral phase transition of QED3. If NJ > 2 then 
the effective low energy theory is chirally broken at zero 
temperature. For the superconductor this means that the 



theory displays a phase transition between the supercon- 
ductingand the antiferromagnetic phase when doping is 
varied [1, Q • If on the other hand Nj < 2 the system 
goes from a superconducting to a pseudogap-phase when 
underdoped. 

The value of Nj has been investigated in a number of 
studies from Dyson-Schwinger eq uations (DSEs) in the 
continuum j,!! [pJ3J|, [IJiaJlll [H, [H [lA] and lattice 
gauge theory [1^, [13, [3 \1m ■ While the former allows 
for numerical as well as analytical studies in principle 
there remains the question of the importance of trunca- 
tion artefacts. For Nj a partial answer has been obtained 
in , where it was found that the details of the fermion- 
photon vertex only have a minor quantitative impact on 
Nj. In particular for all vertex dressings employed Nj 
stayed well above Nf = 2. 

In all studies of dynamical chiral symmetry breaking 
in continuum QED3 a clear separation of scales has been 
found. The intrinsic scale of the theory is given by the 
dimensionful coupling constant a — Nfe'^/8. Then a 
second and much lower scale is given by the dynamically 
generated fermion mass ^(0) in the chiral limit. Related 
to this one has small values for the other order parame- 
ter, the chiral condensate. On the lattice with its finite 
volume this separation of scales is hard to bridge. Conse- 
quently, recent studies for the number of flavours Nf — 2 
and Nf = 4 [l3, determined bounds on the chiral 
condensate, but no definite value for Nj"* could be ex- 
tracted. A definite signal for chiral symmetry breaking 
was obtained only for Nf = 1 [l^ . A very recent calcula- 
tion claims that Nj « 1.5 [l^ with the caveat of volume 
and discretisation artefacts. 

In general, the presence of an infrared cutoff due to the 
finite volume reduces the value of the critical number of 
flavours. This has been demonstrated by Gusynin and 
Reenders in a simple truncation scheme for the DSEs 
(20| . They considered an approximation to the DSE 
for the fermion self-energy which neglects corrections to 
the fermion wave-function and the fermion-photon ver- 
tex. The photon is then given by its leading behaviour 
in the 1/Nf expansion. In this work we elaborate on 
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FIG. 1: The Dyson-Schwinger equations of the photon and 
fermion propagators in diagrammatic notation. Wiggly hnes 
denote photon propagators, straight hnes the fermion. A 
black dot denotes a bare fermion- photon vertex while the open 
circle denotes a dressed one. 



these findings by considering a more sophisticated trun- 
cation scheme which explicitly takes into account non- 
perturbative effects in all of these quantities. In addition 
we follow a different strategy to assess the volume effects 
by evaluating the system on a three-torus^. 

This work is organised as follows. In section |TT] we re- 
call results for QED3 in the infinite volume/continuum 
limit. We present the Dyson-Schwinger equations for the 
fermion and photon propagators and discuss their asymp- 
totic properties. In section IIIII we recall general condi- 
tions for chiral symmetry breaking on a finite volume, 
formulate the DSEs in a box and discuss our numerical 
methods to solve these. In section ITVl we present our nu- 
merical results for the propagators at finite volume and 
the critical number of flavours N'j as a function of the 
box size. In section |V] we conclude. 



II. QED, IN THE INFINITE 
VOLUME/CONTINUUM LIMIT 

A. The Dyson-Schwinger equations in QED3 

We consider QED3 with a four-component spinor rep- 
resentation for the Dirac algebra and Nf fermions. This 
allows a definition of chiral symmetry similar to the cases 
of QED4 and QCD4. With massless fermions, the La- 
grangian has a U{2Nf) "chiral" symmetry, which is bro- 
ken to SU{Nf ) X SU{Nf ) X C/(l) X [/(I) if the fermions 
become massive. The order parameter for this symmetry 
breaking is the chiral condensate which can be deter- 
mined e.g. via the fermion propagator. 

The Dyson-Schwinger equations for the photon and 
fermion propagators in Euclidean space are given dia- 
grammatically in Fig. [TJ They read explicitly 



- Z.Nfe'J ^Tr[7^5(g)r,(g,fc)5(fc)], 



Zie' I S^lt.S{q)T,{q,p)D,,{k), (2) 



(27r)3 



with the momentum routing fc^ = (7^ — . Here Nj de- 
notes the number of fermion flavours and Zi is the renor- 
malisation constant of the fermion-photon vertex T^. 

In Landau gauge the general form of the dressed 
fermion propagator Siji) and the photon propagator 
D^^{p) is given by 



S{p) 



iM{P^)+B{p^) 
p^A^{p^)+B^{p^) ' 



Df^^ip) = (^M"^ 2 ; 2 
* p'^ J p'^ 



= 6. 
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p2 ; p2(^3 + n(p2)) 



(3) 



(4) 



with the photon dressing function G{p^), the photon 
polarisation Ti{p^) and renormalisation function and 
the fermion dressing functions A{p'^) and B{p^). These 
can be rearranged into the renormalisation group invari- 
ant fermion mass function Af(p2) = B{p^)/A{p'^) and 
the fermion wave function Zf{p'^) = 1/A{p'^). Another 
renormalisation group invariant is the 'running coupling' 

G{p^) built from the renormalisation point dependent 
photon dressing function G^p^) and the renormalised 
coupling e^. The bare renormalised fermion propaga- 
tor is given by Sq^{p) = Zi{ij) -I- Zmm) where m is the 
renormalised current fermion mass, Z^n is the mass renor- 
malisation function and Z2 the corresponding one for the 
fermion wave function. Note that in QED we have the 
Ward-Takahashi identity Zi — Z2. 

The gauge dependence of the fermion mass and wave 
function and, correspondingly, the fermion-photon ver- 
tex has been a much debated issue in the past, see e.g. 
[1 m, m, m M, m, m, M, M, hi and references 
therein. There in particular the technical question of how 
to truncate the fermion-photon interaction to obtain a 
gauge covariant pattern of chiral symmetry breaking and 
restoration has been discussed. This issue is no t q uite 
settled yet, however there are clear indications [sl. l32l [33| 
that Landau gauge is preferred in the sense that it allows 
for particularly simple approximation schemes^. This is 
why we prefer Landau gauge in this work. 

A range of ansatze for the fermion-photon vertex 
have been investigated in Q. There it has been found 



1 A corresponding technique has been applied in QCD4 to deter- ^ 1° all other linear covariant gauges parametrised by the gauge 
mine finite volume effects in the quark and gluon propagators parameter ^ a gauge dependent scale appears which compli- 



2II . I22I . I23II cates matters considerably. 
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that the critical number of flavours Nj obtained with 
the most elaborate construction, obeying the Ward- 
Takahashi identity, is almost similar to the one of the 
most simple ansatz, a bare renormalised vertex 

r,. = Zi-f^ . (5) 

Therefore to keep matters as simple as possible we will 
only present results for the bare vertex approximation 
in this work. We did check, however, that more sophis- 
ticated ansatze do not alter the main conclusions pre- 
sented below. In addition, a general analysis of the in- 



frared behaviour of QED3 in the framework of functional 
renormalisation group equations indicates that in Landau 
gauge the bare vertex ([5]) may even be the best possible 
choice [345. 

Substituting this vertex into the fermion and photon 
DSEs, taking appropriate traces and contracting the pho- 
ton DSE with the projector 

v,M = 5,.~(:^ (6) 

with C = 1 we arrive at 



G(p2) 



Zz + n(p2) 



Z, + Zle^ 



Z-i - Zle^Nf J 



(27r)3g2^2(q2)+^2(g2) p 
(2^ ^2^2(^2) _^ ^2(^2) 

(fiq 1 1 



(27r)3 q-^A^{q2) + B^{q^) k^A'^{k^) + B^{k^ 



Wiip'-.q^k^) 



(7) 
(8) 
(9) 



where we used Zi = Z2 and the subtracted kernel 



2P(3-C) 
3p2 



(10) 



with 



Ck^ 



1 - C 2Cg2 



p- 



(11) 



in the photon equation. As explained in appendix A 
of ref. [5i] the subtraction of the term proportional to 
(3 — ^) in Eq. (I10|) is necessary to avoid spurious linear 
divergences in the photon-DSE generated by the regular- 
isation procedure (a hard cutoff) used in the numerical 
treatment of the equations. The choice = 1, i.e. the 
transverse projection of the photon equation, is manda- 
tory to avoid the back-reaction of spurious longitudinal 
terms into the right hand side of the photon equation. 
Compared to 0, where C = 3 has been used, this treat- 
ment leads to a quantitatively improved value for the 
critical number of flavours, see below. Note, however, 
that the effects of varying are quantitatively small in 
general; all qualitative conclusions derived in [5i] and also 
here are independent of the choice of 



B. Asymptotic behaviour of the propagators 

An often used approximation to determine the asymp- 
totic behaviour of the fermion and photon dressing func- 
tions has been the 1/A'^y-expansion. This expansion is 



equivalent to a perturbative expansion for small while 
keeping a = Nje^/S flxed. As QED3 is an asymptotically 
free theory this expansion does provide correct answers 
in the ultraviolet. For the photon polarisation and the 
vector dressing function of the fermion one flnds for Nf 
massless fermion flavours lal 



n(p2 > a) = - 
Aip"^ > a) = 1 



'/ 



a 
P 



(12) 
(13) 



For the mass function of the fermion one can use the 
operator product expansion to obtain 



M (p2 > a) 



2 + e 



(14) 



with the chiral condensate (^'4'). Note that the conden- 
sate can also be determined from the trace of the fermion 
propagator: 



AZn 



d^q 



B{q' 



(27r)3 g2^2(g2) ^_52(g2) 



(15) 



In the infrared momentum region the l/A^/-expansion 
is clearly not sufflcient. Here one has to resort to a self- 
consistent analysis of the DSEs in terms of asymptotic 
power-laws. This method is well developed in QCD4 
In 111,133, Hi] and has been adapted to QEDg in Ref. [1]. 
Such an analysis is valid if no scales are present, i.e. in 
the deep infrared momentum region p <C a and only in 
the absence of fermion masses. The appearance of in- 
frared power laws with potentially fractional anomalous 



4 




truncation scheme. In principle one can reproduce (jl9p 
also in a numerical treatment of the DSEs with a hard 
cutoff when vertex corrections around the cutoff scale are 
included, see for details. However, this procedure is 
much more involved than the simple subtraction scheme 
of Eq. Uni). Since the difference between (HH]) and (HH) is 
rather small in the interesting region above Nf = 2 and 
therefore unimportant for all of the conclusions of the 
present work we resort to the simple subtraction scheme 
PH)) and consequently reproduce (US]) in our numerics. 

The critical number of flavours Nj where chiral sym- 
metry is restored can be determined analytically from 
the DSE for the scalar fermion dressing function B{p^). 
Again we refer the reader for details to Ref. Q and merely 
state the result 



FIG. 2: The infrared exponent k as a function of the number 
of flavours Nf in the symmetric phase of QED3. 



dimensions is thus an indicative and characteristic prop- 
erty of the symmetric phase of QED3 . 

For a bare fcrmion-photon vertex the infrared be- 
haviour of the vector fermion dressing function and the 
photon polarisation in the chirally symmetric phase can 
be written as 



•P 



n(p2) ^ p 



-1-4k 



(16) 



In case of a dressed fermion-photon vertex according to 
the Ward-Takahashi identity the corresponding power- 
laws are 



•P 



n(p2) ^ p 



-1-2k 



(17) 



These expressions solve the DSEs in the chirally symmet- 
ric phase as described in detail in Q (see also [34I for a 
re-derivation from a slightly different perspective). 

The value of k as a function of Nf can then be deter- 
mined from the DSEs using a procedure given in detail 
in The solutions for a bare vertex with and without 
the subtracted term in Eq. (fTO)) is plotted in Fig. [2l For 
Nf > 1.5 the exact solutions can be represented by the 
fits 



sub 
'^bare 



l^ba 



0.142 


0.002 


Nf ^ 


Nj 


0.137 


0.058 


Nf ^ 


' ^1 



0{1/N^f) . 
0{llNj) . 



(18) 
(19) 



A corresponding fit for the bare vertex and a Ward- 
Identity improved vertex for C = 3 is given in The 
difference between (fTS]) and (fT9|l represents the system- 
atic uncertainty in our numerical calculation due to the 
cutoff regularisation procedure used. For dimensional 
regularisation the subtraction of the term in (fTO|) is not 
necessary and therefore the unsubtracted result Eq. (fT51) 
should be viewed as the exact one for the bare vertex 



(A^/)re « 3.84, 

{N'j)bare ~ 3.61. 



(20) 
(21) 



The corresponding result from the functional renormali- 
sation group is {N^)bare ~ 3.6 34] in agreement with our 
result for the unsubtracted equation. Numerical results 
for Ward-identity improved vertices lead to results in the 



range of 3.5—4 [5|. These numbers can be contrasted with 
{Nf)i/Nj = 32/7r2 « 3.24 from the l/7V/-expansion 0. 
Note that all these results are far above Nf = 2 relevant 
for the description of high-Tc cuprate superconductors as 
discussed in the introduction. 

Finally we wish to emphasise that these power law so- 
lutions can only be obtained when the full structure of 
the propagator DSEs is taken into account. The authors 
of Refs. U[l3| did not find a solution corresponding to a 
symmetric phase in their truncation scheme because the 
feedback from the function A onto the vacuum polarisa- 
tion is not taken into account. This then prohibits the 
appearance of power laws and therefore does not allow 
for the appearance of the chirally symmetric phase and 
should be discarded 0, . 



III. QED3 AT FINITE VOLUMES 

A. Chiral symmetry breaking in a box 

Before we embark on our investigation, let us recall the 
gen eral finite volume behaviour of the chiral condensate 
[39j. The fermion propagator in its spectral representa- 
tion is given by 



SAix,y) 



j2 Unjx) uliy) 



m 



(22) 



where Un{x) and A„ are eigenf unctions and eigenvalues 
of the Euclidean Dirac operator, Un{x) = A„M„(a;). 
The gauge field A is treated as an external field. These 
eigenfunctions occur either as zero modes or in pairs of 
opposite eigenvalues. Setting x = y, integrating over x 
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and neglecting the zero mode contributions, one obtains 

^y^-(-^-)---E;;^- (23) 

y A„>0 " 

The chiral condensate can be deduced by averaging the 
left hand side of this equation over all gauge field con- 
figurations and then taking the infinite volume limit to 
give 

oo 

(qq) = -2m / dX (24) 



where p(A) is the mean level density of the spectrum, 
which becomes dense in the infinite volume limit. In the 
chiral limit, m —^ 0, only the infrared part of the spec- 
trum contributes and one finally arrives at the Banks- 
Casher relation 

(qq) = -npiO) . (25) 

If the two limits arc interchanged, i.e. if one takes 
the chiral limit before the infinite volume limit, one 
has a discrete sum in Eq. (|23p and the infrared part 
of the spectrum cannot trigger a non-vanishing chiral 
condensate: chiral symmetry is restored. If, however, 
at a given volume the explicit fermion mass m is not 
too small, one can still observe the spontaneous forma- 
tion of a quark condensate. If the factor (to^ -|- A^)~^ 
varies only slightly with n, the sum in Eq. (|23p can still 
be replaced by an integral and Eq. remains valid. 
For this to be a legitimate approximation one needs 
771 > AA l/Vp{X) = n/{V\{qq)\), at the lower end 
of the spectrum. Thus one obtains the condition [s^ 

Vm\{qq)\:^TT. (26) 

This relation reveals the crux of the matter. In principle, 
if the volume is large enough, the necessary quark masses 
are academically small and may even be neglected in the 
numerical treatment. What counts as large in this con- 
text, however, depends sensitively on the size of the chiral 
condensate. In QED3 the condensate becomes extremely 
small already well below the critical number of flavours 
of the chiral phase transition Q. Thus although a 
formulation on a finite volume may do well for Nf = 1, 
all signals of dynamical chiral symmetry breaking will be 
lost already well below the of the theory in the in- 
finite volume/continuum limit. This behaviour will be 
quantified below. 

B. The DSEs on a torus 

On a compact manifold, the photon and fermion fields 
have to obey appropriate boundary conditions in the time 
direction. These have to be periodic for the photon fields 
and antiperiodic for the fermions. For computational rea- 
sons it is highly advantageous, though not necessary, to 




FIG. 3: Two-dimensional sketch of the momentum grid dual 
to the three-torus for a fixed Cartesian momentum cutofi'. 
The hyperspheres depicted by dashed lines are not complete in 
the sense that additional momentum points on these spheres 
are generated if the cutofi' is increased. The torus equiva- 
lent of an 0(3)-invariant cutofi^ used in our calculations sums 
only over complete hyperspheres, which are indicated by fully 
drawn circles. 

choose the same conditions in the spatial directions. We 
choose the box to be of equal length in all directions, 
Li = L2 = Lj, = L, and denote the corresponding vol- 
ume V = . Together with the boundary conditions 
this leads to discretised momenta in momentum space. 
Thus all momentum integrals appearing in the Dyson- 
Schwinger equations are replaced by sums over Matsub- 
ara modes. 

On a torus with antiperiodic boundary conditions for 
the fermion fields, the momentum integral changes into 
a sum of Matsubara modes, 

^ E O'O. (27) 

counting momenta q„ — 3(27''/i)(?T.i -I- l/2)ei, 

where ei are Cartesian unit vectors in Euclidean mo- 
mentum space. For the photon with periodic bound- 
ary conditions the momentum counting goes like qn = 
3(27r/i)(ni)ei. For the numerical treatment of the 
equations it is convenient to rearrange these summations 
such that they represent a spherical coordinate system 
[2l| . see Fig. [3] for an illustration. We then write 

E (•■■) -isE^--)' (28) 

where j counts spheres with qnqn — const, and m num- 
bers the grid points on a given sphere. The corresponding 
momentum vectors are denoted (\m,j and their absolute 
values are given by j = llm.j I • It is then a simple mat- 
ter to introduce the torus equivalent of an 0(3)-invariant 
cut-off by restricting j to an interval [1,N]. This proce- 



6 



dure is equivalent to 'cutting the edges' of our torus as 
indicated in Fig. [31 



b{pI) 



1 



Note that the momentum argument k of the photon self- 
energy in the fcrmion DSE is the difference ki^ij^q — 
Pi,i ~ Qj.q of two antiperiodic Matsubara momenta and 
thus lives on a momentum grid corresponding to peri- 
odic boundary conditions, as it should. 

The DSEs can be solved numerically employing well es- 
tablished methods. Our numerical method on the torus 
is outlined in Ref. (2^ . the corresponding continuum 
method as well as details on the renormalisation pro- 
cedure of the DSEs are given in Ref. [§]. 

Note that the propagators determined from the contin- 
uum version of the DSEs, Eqs. are independent of 
the regularisation procedure. In our numerical calcula- 
tions in the infinite volume/continuum limit we use a sub- 
tracted version of these equations and an 0(3)-invariant 
UV-cutoff which can be sent to infinity at the end of each 
calculation. These DSEs therefore represent not only the 
infinite volume limit but also the continuum limit (in co- 
ordinate space) of the representation given by Eqs. ([THS]) 
of the DSEs on a torus. We use the phrase infinite vol- 
ume/continuum limit to indicate this simultaneous re- 
moval of both an ultraviolet and an infrared cutoff. 



IV. NUMERICAL RESULTS 

A. Finite size effects: fixing the UV cutoff' 

Before we investigate the finite volume effects of vari- 
ous quantities of interest we have to clarify whether there 
are sizeable effects due to the finite size of the system cor- 
responding to a fixed ultraviolet momentum cutoff. Part 
of these effects are removed by our choice of 'cutting the 
edges' of the torus, described in the previous section. The 
remaining finite size effects can be evaluated by varying 
the size of the cutoff. Here, a natural minimal cutoff is 
given by the intrinsic scale a = Nf e^/8 of QED3. Thus, 
working in the range < Nf < 8, we anticipate that fi- 
nite size effects are of minor importance for cutoffs larger 



The resulting DSEs are then given by 



I 

than the intrinsic scale, i.e. > e^. Indeed, this is the 
case as can be seen from the two plots in Fig. 21 Shown 
are the fermion mass function MijP') and wave function 
Zf{p'^) as a function of scaled momentum /e^ for var- 
ious cutoffs between A^y = 0.15 e** and A^y = 39.5 e**. 
The number of fermion flavours is chosen to be Nf = 1, 
similar results are obtained for other choices. The box 
volume for the fermions is related to the lowest momen- 
tum point by y = L'^ with L = \/iTT /pmin- Here we 
use L = 390/e^. Note that for the photon the box vol- 
ume is related to the lowest momentum point \yy V = 
with L = 27r/p„i„. The same box length then results in 
a slightly different value for the lowest momentum point 
in the photon dressing function. We refrain from showing 
the photon explicitly here, since the finite size effects are 
similar to the ones for the fermions. 

In Fig. m the finite size effects are visible only for the 
smallest cutoff. There are effects in both the infrared and 
ultraviolet momentum regions, where a number of mo- 
mentum points deviate from the results with larger cut- 
offs. This is true for the renormalisation point dependent 
fermion wave function, normalised such that they match 
the continuum results, but also for the renormalisation 
point independent fcrmion mass function. The perhaps 
surprising observation that the variation of an ultraviolet 
cutoff also affects the infrared behaviour of the dressing 
functions implies a certain entanglement between the in- 
frared and ultraviolet modes of a gauge theory. Similar 
effects have been found in four dimensional Yang-Mills 
theory, see e.g. [s^ and references therein. 

We also need to comment on the fact that we observe 
dynamical chiral symmetry breaking in our system at all, 
despite in practice working with a vanishing bare quark 
mass. In the continuum formulation this entails working 
in the chiral limit. However, due to the formal reasons 
discussed in section UlI Al this cannot be true on a finite 
volume, since taking the chiral limit before the infinite 
volume limit inevitably leads to the loss of dynamical 
chiral symmetry breaking. However, this is not what 
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FIG. 4: Fermion mass function M{p^) and wave function 
Zflp^) as a function of scaled momentum p^/e*. The volume 
of the box corresponding to the smallest available momen- 
tum is kept fixed and the ultraviolet cutofi^ is varied. Here 
counts the number of momentum points in one direction 
of our torus, i.e. N = 100 means that we employed a torus 
with 200^ momentum vectors. This translates into the cutoffs 
Afjy = 0.15 e'', 0.62 e*, 2.47 e*, 9.86 6*, 39.5 e^. 

we observe here. Indeed, the volumes we use are large 
enough to allow for extremely small bare fermion masses 
according to the relation of Eq. (|^. In practice, this 
allows one to neglect the fermion mass entirely in the 
numerical treatment of the DSEs. However it is impor- 
tant to keep in mind that this means we are working close 
to the chiral limit, but not in the chiral limit. 

We conclude this subsection with the observation that 
from about A^y = onwards almost no finite size ef- 
fects are present. We therefore consider this value a lower 
bound for admissible cutoffs on a torus, in agreement 
with our general considerations concerning the natural 
scale in QED3. This finding also agrees with correspond- 
ing results on finite size effects in lattice simulations, see 
and Refs. therein. In order to be absolutely sure 
that cutoff effects do not play any role in what follows 
we use the somewhat larger cutoff Afjy — 2.35 from 





10"' 10"' 10"^ 10"' 10"^ 10"' 10° 10' 
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FIG. 5: Fermion mass function M{p^), wave function Zf{p^) 
and photon self energy n(p^) as a function of scaled mo- 
mentum p^/e*. The ultraviolet momentum cutoff is kept 
fixed at AfjY = 2.35 and the box length is varied from 
Le^ = 200 - 3000. 



now on. 
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B. Finite volume effects: towards the infinite 
volume limit 

We are now in a position to study the finite volume 
effects occurring for the fermion and photon propagators 
on a torus. To this end we keep the ultraviolet cutoff of 
our system fixed and vary the infrared cutoff in the range 
of L = 200 — 3000. We show the resulting behaviour of 
the fermion mass and wave function as well as the pho- 
ton together with the corresponding continuum results in 
Fig. [5l Again, we choose Nj — 1. The variation of the 
volume clearly results in the loss of a substantial amount 
of generated fermion mass when the volume gets smaller 
and smaller. For even smaller volumes than shown in 
the figure chiral symmetry is restored in agreement with 
the condition Vm\{qq)\ 3> tt discussed in the subsections 
IIII Al and IIV Al On the other hand, we observe that ex- 
tremely large volumes are needed to account for the full 
effect of dynamical chiral symmetry breaking observed in 
the infinite volume/continuum limit. This is in marked 
contrast with the behaviour of the quark sector of QCD4 
[m and can be explained by a closer look at the scales 
involved. Whereas in QCD4 the generated quark masses 
{M (0) 300 - 400 MeV) are of the same order as the 
intrinsic scale of the system {-^-qcd ~ ^50 MeV) , the 
situation is clearly different in QED3. Here our char- 
acteristic scale is of the order a = whereas the 
generated fermion masses are of order 10~^ e^, as can be 
seen from the plot. To keep volume effects small, this 
scale has to be well accommodated by the system on a 
box, which translates to a lowest momentum to be much 
smaller than p'^ = 10~^ e^. Indeed, this is what we ob- 
serve. Choosing the volume large enough that the lowest 
momentum is well below this scale we approach the in- 
finite volume/continuum limit. This is true for all three 
dressing functions, the fermion mass function M(p^), the 
fermion wave function Z f(p-) and the photon self energy 

n(p2). 




FIG. 6: Upper plot: Fermion wave function Zf{p^) as a func- 
tion of scaled momentum p^/e* and the number Nf of fermion 
flavours. Note that for clarity of the figure 

frompVe" = 10"^ 

upwards only results for a selected number of momentum 
points on the torus are shown. Lower plot: The resulting 
value of the infrared exponent function of Nf for dif- 

ferent volumes. 



C. Finite volume effects: Zf{Nf) and K,{Nf) 

Let us now examine the infiuence of finite volume ef- 
fects on the critical value Nj, where the system undergoes 
a phase transition from the chirally broken into the chi- 
rally symmetric phase. This phase transition is marked 
by the change of the infrared asymptotics of the fermion 
wave function for Nf < N'j this function is a con- 

stant in the infrared, whereas for Nf > N'ijt develops a 
power law with iV/-dependent exponent k ■ 5!| . In a sense, 
K can be viewed as an order parameter of this phase tran- 
sition"^ . We exhibit this behaviour of k by fitting a power 



^ Note, however, that the phase transition is not a first or sec- 
ond order transition but has similar properties to the conformal 
transition of QED in four dimensions [3, |43I . 



law in the infrared to the fermion wave function 

Zf{p^)^C{p^)-\ (32) 

which is related to the fermion vector dressing function 
by Zf = \/A. Both, the coeflacient C and the power k 
may depend on Nf. 

To visualize this procedure the fermion wave func- 
tion Z is shown in the upper panel of Fig. [5] for a 
box with Le^ ~ 800 for different numbers of flavours 
Nf. One clearly sees the aforementioned behaviour: for 
Nf < 2.0 the function develops a constant in the in- 
frared region {i.e. k — 0), whose value is proportional 
to Nf. The phase transition occurs for this specific vol- 
ume for 2.0 < < 2.2. For Nf > 2.2 one observes a 
power law of the wave function at small momenta. The 
coefficient of this power law is still proportional to Nf, 
however the exponent k decreases with growing Nf ac- 
cording to Eq. I|18p. Consequently we observe a decrease 
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of the function Zf{p'^) with growing Nf in the symmet- 
ric phase. As a result one could determine Nj as the 
number of flavours for which Zf is maximal at a given 
infrared momentum. This is also the case in the infinite 
volume/continuum limit [5]. 

As an (equivalent) alternative we determine Nj by 
fitting the power law (|32p to our numerical results for 
Zf{p^) in the infrared momentum region. There is a 
caveat here: similar to four-dimensional Yang-Mills the- 
ory one observes that the power law ((32)) can only be seen 
for momenta 1/L ^ p <C e^. This behaviour is generic 
on a torus [1^. To obtain significant results for our in- 
frared coefficients we therefore have to perform the fit 
for momentum points significantly larger than the lowest 
one. In practice we chose momenta from the fourth point 
onwards. 

The resulting dependence of the exponent k on Nf is 
plotted in the lower panel of Fig. [5] for several volumes 
of the box together with the corresponding function in 
the infinite volume/continuum limit. At infinite volume 
the function K{Nf) is zero for Nf < N'j k, 3.84 and equal 
to the analytic result of Eq. ([^D]) in the symmetric phase 
above Nj. For a given finite volume the phase transition 
is still indicated by a maximum in k,{N f). However, there 
is an additional region at Nf < Nj{L), where the func- 
tion rises slowly towards its maximal value. This region 
is generated by the finite infrared cutoff of the system in 
a box, which prevents the fermion wave function Zflp^) 
from bending towards a finite asymptotic value at = 0. 
This effect mimics a power law at values of Nf where the 
system is still in the chirally broken phase. With increas- 
ing volume this region gets smaller and smaller until it 
reaches the sharp rise of the function in the infinite vol- 
ume/continuum limit as shown in Fig.[6l 

As a result we find a critical number of flavours which 
depends upon the volume of the torus. The explicit val- 
ues are shown in Fig. [7] together with corresponding re- 
sults from Ref. [2^. In the following we concentrate on 
the solutions close to the chiral limit and postpone the 
discussion of the (orange) curve with large bare fermion 
mass to the next subsection. 

For volumes that are not too large, our solutions and 
the ones reported in Ref. [13] can be well fitted by a form 



N^f ^a-b/{Le^) 



2x1/3 



(33) 



For our results we obtain a — 3.23 and b — 10.64, shown 
as dashed line in the plot. Interestingly, to high preci- 
sion the exponent of the box length in this fit is given 
by —1/3, although we do not have a good explanation 
for exactly this behaviour. In terms of (1/L)^/'^ the fit 
then suggests a linear extrapolation to the infinite vol- 
ume limit. However, it turns out that this linear be- 
haviour breaks down for extremely large volumes. This 
is evident for the results from Gusynin and Reenders [20[ . 
The deviation from the linear behaviour is quantitatively 
important: for the truncation of Gusynin and Reenders 
a linear extrapolation to the infinite volume results in 



Ref. [20] ' 
linear extrapolation 
'chiral' 

linear extrapolatioiL 
m = 0.0024 e' 




8000 



1/3 



1000 



1/3 



300 

2-1/3 

(Le ) 



1/3 



125 



1/3 



64 



■1/3 



FIG. 7; The critical number of flavours A'^^ as a function of 
inverse box length In the 'chiral' limit we show the 

results of our calculation compared to the corresponding one 
of Ref. [20| ] in a simpler truncation scheme. The dashed lines 
are linear extrapolations to the infinite volume limit. We also 
show results for a bare fermion mass comparable to the ones 
used in lattice calculations. 



N 



f 



N 



2.52 as compared to the infinite volume result 
y — 3.2. (Note that the curve gets extremely steep 
for the very largest volumes, which are not plotted.) For 
our results we find a value of Nj « 3.23 for the linear ex- 
trapolation, whereas the analytical infinite volume result 
is N^ « 3.84. Though this difference is not huge, it high- 
lights the need for non-linear extrapolation procedures to 
the infinite volume limit. Unfortunately our largest vol- 
umes are not yet large enough to penetrate this region of 
nonlinearity. 



D. Non- vanishing bare fermion masses 

In this subsection we investigate the consequences of 
explicit bare fermion masses toq in the Lagrangian of 
QED3 on the chiral phase transition. To this end we 
first determined the properties of the system in the infi- 
nite volume/continuum limit by solving the correspond- 
ing DSEs for a range of different bare fermion masses. 
For any given mass and varying Nf we find that the 
fermion wave function never develops a pure power law 
in the infrared. Instead there is a region < p < A/^, 
where the function is constant and a region Am < p < a 
where a power law is present. The scale Am depends 
on the explicit fermion mass as well as on the number 
of flavours. We conclude from this that away from the 
chiral limit the infrared exponent k is no longer an order 
parameter. As concerns the global behaviour of Zf{p'^) 
we still observe the behaviour discussed above Eq. (|33p : 
at a small enough momentum p^ the value of the function 
Z f{p'^) increases with N f up to a certain point at a criti- 
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FIG. 8: The chiral condensate as a function of the fermion 
mass for a range of volumes at Nf = 1.5. SmaU volumes give 
rise to a vanishing chiral condensate and are nearly degener- 
ate, falsely indicating weak finite volume effects. 

cal value Nj and decreases again if Nf grows further. In 
the chiral limit this critical value marked the chiral 
phase transition. Here, however, this seems not to be so. 
The fermion mass function il/(0) decreases exponentially 
with Nf without a trace of a rapid change around Nj. 
We therefore conclude that there is no phase transition 
in the infinite volume/continuum limit for QED3 with 
non- vanishing bare fermion masses ttiq. 

On a torus with a given volume, however, the scale 
A7fl(7V/) decreases with Nf and can become lower than 
the lowest momentum point Pmin on the torus. Conse- 
quently one then sees a pure power law in Zf{p'^) for 
momenta Pmin < p < a. The critical number of fermion 
flavours Nj where this transition is observed coincides 
with the value of Nf where Zf{p'^^^) is maximal. Thus 
in a sense one observes a chiral transition on a torus even 
for non-vanishing toq when there is none in the infinite 
volume continuum limit. For the choice m = 0.0024 
the resulting values of iV^ are plotted against (ie^)^/"^ in 
Fig- El Apart from our smallest volumes the resulting val- 
ues of Nf are significantly larger than in the 'chiral' limit. 
We also observe that in the infinite volume limit the crit- 
ical number of flavours goes to infinity, in agreement with 
our findings discussed in the previous paragraph. Away 
from the chiral limit QED3 on a torus shows a chiral 
phase transition when there is none in the inflnite vol- 
ume/continuum limit. This implies that extreme care is 
needed when one attempts to extract information on Nj 
from lattice calculations with finite fermion masses. 

Nevertheless it is interesting to compare to lattice re- 
sults. As explained above the important scale in assessing 
finite volume effects is the lowest momentum point avail- 
able on the torus compared with the generated fermion 
mass. In our setup with antiperiodic boundary condi- 
tions in space and time directions this scale is given by 
Lose — ^/?>t^ /Pmin- On the lattice one usually imple- 
ments antiperiodic boundary conditions in the temporal 



direction and periodic ones in the two spatial directions. 
This results in Liatt = T^/Pmin on the lattice. The cou- 
plings are related by ~ l/(/?a), where (3 is the dimen- 
sionless inverse coupling on the lattice and a the lattice 
spacing. With ~ 1 and (3 ~ C(l) one then obtains 
Ldses^ — VSLiattP- This means that we should com- 
pare the results of contemporary lattice calculations on 
80'^-lattices with our values for L « 140. From the plot of 
Fig. [7] we then find chiral symmetry breaking for Nf = 1, 
whereas at Nf = 1.5 the system is in the chirally symmet- 
ric phase. This is in agreement with recent lattice results 
[HI and shows that these results are compatible with our 
value of Nf « 3.8 in the infinite volume/continuum limit. 

Finally we investigate the behaviour of the chiral con- 
densate on a torus as a function of the explicit fermion 
mass. To this end we extract the condensate from our 
fermion mass function with the help of Eq. (jlSp at a 
given fixed ultraviolet momentum cutoff. Our results 
for Nf = 1.5 and different box lengths L are shown in 
Fig. [51 Clearly, for small volumes the condensate de- 
creases linearly with decreasing bare fermion mass toq 
and extrapolates to zero in the chiral limit. For large 
enough boxes this behaviour changes and we find a fi- 
nite value in this limit in agreement with results from 
our continuum DSEs. It is interesting to note that the 
curves for small volumes are almost degenerate. A vol- 
ume analysis in this region would therefore indicate weak 
finite volume effects where in fact there are large effects 
when the volume is increase further. Again, this result 
should serve as a caveat for the interpretation of lattice 
results. 



V. SUMMARY 

In this work we investigated volume effects on the chi- 
ral phase transition of QED3 as a function of the number 
of flavours, Nf. To this end we solved a coupled system 
of Dyson-Schwinger equations for the fermion and pho- 
ton dressing functions in the inflnite volume/continuum 
limit and on a three-torus. We worked in a truncation 
scheme that in the infinite volume/continuum limit re- 
produces a critical number of flavours Nf « 3.61 — 3.84, 
a number close to the one obtained with more involved 
approximation schemes 

Examining the same system on a torus we found con- 
siderable volume effects. These can be explained due 
to the presence of scales of vastly different magnitude 
in QED3. On the one hand one has the natural scale 
a/e'^ = Nf/8, which is of order one. On the other hand, 
the dynamically generated fermion masses are orders of 
magnitude smaller. For example, one has M(0) = 0.05e^, 
for Nf = 1 as can be seen from Fig. [5l While the first 
scale remains of the same order, the second one rapidly 
decreases when the number of flavours becomes larger. 
When the generated fermion mass drops below the in- 
frared cutoff of the system, given by the inverse of the 
box length, chiral symmetry breaking disappears and the 
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system falls into the chirally symmetric phase. We quan- 
tified these effects and presented results for the critical 
number of fiavours as a function of the box length. They 
agree qualitatively with corresponding results of [20] , de- 
termined in a simpler truncation scheme and a differ- 
ent set-up (continuum limit with finite infrared cutoff). 
Quantitative differences are small. 

Our results confirm the notion that lattice calcula- 
tions at Nf = 1.5 or Nf = 2 need very large physical 
volumes to see dynamical chiral symmetry breaking, let 
alone quantify their finite volume effects. Contemporary 
lattice results cannot yet accommodate for these and con- 
sequently find a system in the chirally symmetric phase 



in quantitative agreement with our findings. We have 
shown that these results are nicely compatible with a 
critical number of flavours w 3.61 — 3.84 in the infi- 
nite volume/continuum limit. 
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